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The vortex hopping motion in a type-II superconductor determines the current-carrying ability 
and consequently the application fields of the superconductor. However, it is not clear how the 
vortices hop between the different pinning regions in the superconductor. Here we proposed that 
there should be magnetic flux exchange between two contacting pinning regions. A system of 
differential equations was constructed to described the flux exchange phenomenon. The numerical 
solutions and approximate solutions of the system were obtained. The results shows that the flux 
exchange reduces the internal field in the weak pinning regions, but increases the internal field in 
the strong pinning regions. Moreover, the flux exchange phenomenon is strongly influenced by the 
superconductor's geometrical size. 



PACS numbers: 74.25.Wx,74.25.0p 

I. INTRODUCTION 

The vortices in a type-II superconductor are gener- 
ally pinned down by pinning centers But there is a 
probability for the vortices to spontaneously hop between 
the adjacent pinning centers due to the vortex thermal 
fluctuation?"'" The hopping motion can be described by 
the Arrhenius equation, which shows that the hopping 
frequency is strongly related to the activation energy of 
the vortices.^ The hopping motion causes the magnetic 
flux (or current) in the superconductor to reduce with 
increasing time, which is usually referred as flux relax- 
ation. This phenomenon exhibits various time evolution 
behaviors because of the complicated field (or current) 
dependence of the activation energyPHSJ Due to the dis- 
sipation of the hopping motion, it determines the super- 
conductor's curr ent carry ing ability and consequently its 
application fields!^ 16 * 17 ! Thus, a larg e num ber of exper- 
imental and theoretical research work a 1 * 3 * 16 " 24 ! have been 
carried out to study the vortex hopping motion. 

In the present work, we wish to study the mutual 
vortex hopping motion between the contacting regions 
with different pinning ability, which we shall call flux ex- 
change. More specifically, a real superconductor may be 
divided into smaller regions whose sizes are larger than 
the superconductor's penetration depth. Because of the 
different pinning ability of these smaller regionsp^l the 
vortices close to the interface between two contacting re- 
gions will hop from one region into another regiorPSHSS 
and vice versa. In other words, the two contacting pin- 
ning regions exchange vortices (flux) . Because of the dif- 
ferent pinning ability, the vortex hopping motions in each 
smaller regions need to be described by different equa- 
tions. This indicates that a flux exchange process is de- 
termined by the parameters of both contacting pinning 
regions. It can not be simply described by the conven- 
tional flux relaxation theory. Thus, we need to consider 
the flux exchange seriously and construct new mathemat- 
ical equations to describe it. 

However, the flux exchange is usually obscured by the 
dominant flux relaxation and is then ignored in magneti- 



zation (or intern al field) measurements carried out over 
the entire sample! 1 ^ 29 * 30 ! To study the flux exchange phe- 
nomenon, we may consider the fact that the vortices in 
a weak pinning region have a higher energy and higher 
hopping frequency, but the vortices in a strong pinning 
region have a lower energy and lower hopping frequency. 
Consequently, there should be more vortices hop from 
the weak pinning region into the strong pinning region 
in a unit time interval. It will result in a nonzero vor- 
tex migration from the weak pinning region to the strong 
pinning region, which changes the vortices quantity and 
consequently the flux relaxation law in each pinning re- 
gion. This indicates that it is possible to study the flux 
exchange by investigating its influences on the flux relax- 
ation behavior. 

In this work, the physical principle of the flux exchange 
phenomenon was analyzed. First, we constructed a sys- 
tem of differential equations to describe the flux exchange 
phenomenon by investigating its influence on the flux re- 
laxation behavior. Next, we calculated the numerical 
results of the system to view the influence of the flux ex- 
change on the flux relaxation. Finally, we calculated the 
analytical solutions of the system under a special approx- 
imation. The flux exchange between the bulk pinning 
region and surface pinning region of a type-II supercon- 
ductor was also discussed. 



II. DIFFERENTIAL EQUATIONS OF FLUX 
EXCHANGE 

As mentioned before, the flux exchange here means the 
mutual vortex hopping motion between different pinning 
regions. In this section, we shall first consider a flux 
exchange process that only involves two pinning regions 
(with sizes larger than the superconductor's penetration 
depth): the p-th pinning regions and q-th pinning regions 
(see Fig. 1). Later, we shall generalize it into a flux ex- 
change process that involves N different pinning regions. 

To describe the flux exchange phenomenon, we need to 
construct a system of differential equations. This can be 
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FIG. 1: (Color online) Schematic diagram of magnetic flux 
exchange between the p-th pinning region and q-th pinning 
region of a type-II superconductor. The black disks represent 
vortices. The arrows represent the hopping directions of the 
vortices. A p is the cross area perpendicular to the magnetic 
field in the p-th pinning region. A q is the cross area perpen- 
dicular to the magnetic field in the q-th pinning region. L pq 
is the line length on the interface perpendicular to magnetic 
field. 



achieved by considering the rate of change of the mag- 
netic flux in each pinning region , which is proportional to 
vortex hopping frequency! 7 * 14 ! 31 ^ The hopping process can 
be described by the Arrhenius equation v = h>Qe~ u l k7 ', 
where v is hopping frequency, vq is attempt frequency, U 
is vortex activation energy, k is Boltzmann constant, and 
T is temperature. Multiplying v by the aver age hopping 
distance s, we obtain the vortex velocity^ 3 ^ 

v = sv e- U/kT . (1) 

In the p-th pinning region, the magnetic flux variation 
in a time interval At includes three parts: (a) Due to flux 
exchange, the vortices in the p-th pinning regio n (cl ose 
to the interface) hop into the q-th pinning regiorPEU by 
— A$ p9 = — 1 /2B p L pq VpAt, where B p is average internal 
field (vortex density) in the p-th pinning region, L pq is the 
line length on the interface perpendicular to field, Vp is 
the vortex velocity in the p-th pinning region (see Eq. ([T]) ) . 
The number j% accounts for the random hopping motion 



in both directions, (b) Due to flux exchange, the vortices 
in the g-th pinning region hop back to the p-th pinning 
region by A$ 9p = 1 / 2 B q L qp v q At, where B q is average in- 
ternal field in the q-th pinning region, L qp — L pq , and v q 
is the vortex velocity in the q-th pinning region, (c) Due 
to flux relaxation (dissipative), the vortices in the p-th 
pinning region decrease by — A$ pr = — e p A p e~ u p/ kT At, 
where e p is a proportional constaniPS and A p is the cross 
area perpendicular to the magnetic field in the p-th pin- 
ning region. 

The total change of magnetic flux in the p-th pinning 
region is now A$ p = — A<5> pr — A$ pq + A<5> 9P . Dividing 
both sides by A p At, we obtain a differential equation for 
the average internal field in the p-th pinning region, that 
is, 

^ - -e e - u v/ kT - e B (2) 

where e pq — 1 /2L pq svo/A p and 

Bpq = B p e- u ^ kT - B q c~ u «/ kT . (3) 

is a function related the flux exchange (Eq.Q is used). 

In the g-th pinning region, similarly, the total change of 
magnetic flux is A$ 9 = — A$ 9r — A$ 9P + A$ pg , where 
— A$ gr = — e q A q e~ u i/ kT At, e q is a proportional con- 
stant, and A q is the cross area perpendicular to the mag- 
netic field in the g-th pinning region. Dividing both sides 
by A q At, we obtain a differential equation for the average 
internal field in the g-th pinning region, that is, 

- _ e e -ujkT _ B (4) 

where e qp = l / 2 L qp si'( ) /A q and B qp = -B pq . 

Eq. ([2]) and Eq. Q constitute a nonlinear system of dif- 
ferential equations that describes the flux exchange be- 
tween the p-th pinning region and g-th pinning region. 
If the flux exchange process includes N different pinning 
regions, then the system of differential equations should 
be generalized into 
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For simplicity, we can simply write system ^ as: 



dBj 
dt 



-Ui/kT 



N 



where e, is a proportional constant and = 
1 / 2 L ij su /A i , Bij = B ie - u */ kT - B je - u i/ kT , and i = 
1, • • • ,JV. The ' on j means j i, i.e., the summation 
does not include en terms. If the i-th pinning region and 
j-th pinning region are not contacted, then their interface 



length Li 



0. Consequently, ej 



and 



the corresponding exchange terms vanish automatically 
(Note that generally ^ eji unless Ai = Aj). 

From mathematics we know that there is no general 
method available for calculating the exact solutions of 
a nonlinear system.^ Thus, we shall analyze the non- 
linear system using numerical methods and approximate 
methods. For simplicity, we shall discuss a flux exchange 
process that only involves two pinning regions Eq.([2]), 
Eq.@. 



III. NUMERICAL SOLUTION 

To carry out a numerical calculation of Eq.([2j), Eq.Q, 
we need the initial conditions B p (0) = B p q, Bq~{0) = B q o 
respectively. Let us now rewrite the system as 



1 -pq 



dB P - , p -U p /kT 
dt ~ €P 

dB q = 
dt 

B p (0) =Bpo,B q (0) =B q0 



-U q /kT , ( 
T t qp y 



B„ 
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-U p /kT 



B q e 



B q e 



-UJkT 



-U q /kT 



(6) 

Furthermore, we still need the detailed expressions of 
the activation energies U p and U q . A number of activa- 
tion energies were proposed in early studies pHIH in prin- 
ciple, we can choose any one of these activation energies 
for both U p and U q . If the number of activation energies 
is N, then the number of systems of differential equations 
is iV 2 . This is a big number and it is impossible to ex- 
haust all the possible combinations in the present paper. 
To obtain a direct view on the effects of flux exchange, 
we shall do numerical calculations using the following two 
combinations: 

(a) . Linear activation energy^ U p (B p ) 
B p /B p o), U q (B q ) = U q0 (l — B q /B q0 ). 

(b) . Logarithmic activation energy^ 
U p0 ln(B pa /B p ), U q {B q ) = U q0 ln(B q0 /B q ). 

The parameters were put as follows: B 



= U p0 (l - 
U p (B p ) = 



po 



B, 



„o 



20 

G, U p o/kT = 15, U q „/kT = 75, e p = e q = 0.6 G/s, and 
e qp —e pq —0, 0.05, and 0.50 s" 1 respectively. In the nu- 
merical calculations, we also noticed that the step of t 
must be very small. Otherwise, the results are unsta- 
ble. It indicates that system ^ is a stiff system. The 
calculated results are shown in Fig 2. 



Fig. 2 shows that the flux exchange reduces B p (t) (the 
internal field in p-th pinning region) , but increases B q (t) 
(the internal field in the q-th pinning region). It looks 
that the g-th pinning region attracts vortices from the 
p-th pinning region. In fact, there is no attracting force 
in the q-th pinning region. It is just because we have as- 
sumed that the activation energy U p < U q , and therefore, 
the vortex hopping frequency v v > v q . This difference re- 
sults in a "deficit" of flux exchange in the p-th pinning 
region, but a "surplus" of flux exchange in the g-th pin- 
ning region. 
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FIG. 2: (Color online) Numerical results of magnetic flux 
exchange. In the calculation, we used the following two ac- 
tivation energy combinations for the p-th and g-th pinning 
regions: (a) Linear activation energies U P (B P ) = U p o{l — 
Bp/Bpo) and U q (B q ) = U q o(l — B q /B q0 ). (b) Logarithmic 
activation energy U P (B P ) — U p o \n(B p o/B p ) and U q (B q ) = 
U q o\n(B q o/B q ). The parameters were put as: B p o = B q o = 
20 G, U p0 /kT = 15, U q o/kT = 75, e p = e q = 0.6 G/s, and 
£ pq =£ qp = 0, 0.05, and 0.50 s _1 respectively. 



IV. APPROXIMATE ANALYTICAL 
SOLUTIONS 

In this section, we shall explore the analytical solutions 
of Eq.([2]), Eq.Q under the following approximation: the 
effect of flux exchange is ignorable for the average inter- 
nal field in the p-th pinning region B p (t), but is significant 
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for the average internal field in the q-th pinning region 
B q (t). Using mathematical language, we can write out 
the approximation as: (a) e pq B pq ss 0; (b) \e qp B qp \ >> 0. 

Later, we shall show that the above approximation can 
be used to describe the flux exchange process between 
the bulk pinning region and surface pinning region of a 
type-II superconductor. 



A. Solutions in the p-th pinning region 



(Note that B pq ^ 0. Otherwise, B q 



B pq — 



and the second approximation (b) \e qp B qp \ » cannot 
be satisfied) . This occurs when the area of the p-th pin- 
ning region A p is so large that e pq — 1 / 2L pq svQ/ A p s» 0. 
Therefore, Eq.([2| reduces to 



dBp 
dt 



-U p /kT 



(7) 



Eq.Q is a decoupled differential equation. It has ex- 
act solutions with the following two activation energies^ 
(with initial condition B p (0) = B p o): 
(1). Linear activation energy: U p (B p ) — U p q(1 — 
Bp/Bpo). The solution is 



B p (t) — B p v 



(8) 



where /3 p = U pQ /kT and r p = B p0 /(f3 p e p ). 
(2). Logarithmic activation energy: 
Upo In (B P Q /Bp), The solution is 



U p (B p ) 



B p (t) — B p0 I 1 



where 7 = l/(/3 p — 1) and t p = B p oj/e p . 



(9) 



1. B qp « B q e- u «' hT 
Substituting B qp sa B q e~ u ?/ kT into Eq.Q, we have 

^ = -e g e- u «/ kT -e q pB g e- u «/ kT . (11) 



Substituting Eq.(10) into Eq.(ll), we have 



e x dx 
x dt 



The first approximation, (a) e pq B pq sa 0, indicates that where 



x = P q B q + (T q e qp ) \ 

X = X(0) = /3 g + (TgCqp)" 1 , 

p q = U q0 /kT, (3 q = P q /B q0 , and r q = l/0 q e q ). 



(12) 



(13) 



(14) 



Integrating both sides of Eq.(12| with respect to t, we 
have 



ln\x\ + J2 

n=l 

where the formula 



(-*)" 



n ■ n\ 



(egpe-^t + d (15) 



dx = In I 



E 



(- x y 



is used and 



d^ln\x \ + Y, 



(-*o) n 



is an integrating constant that is determined by the initial 
condition x(0) = f3 q + (t^)" 1 (or B q (0) = B q0 ). 

The internal field in the g-th pinning region, B q (t), is 
then determined by Eq.((l3|, Eq.(Jl4|), and Eq.(fl5}. 



B. Solutions in the q-th pinning region 

The second approximation, (b) \e qp B qp \ >> 0, indi- 
cates that \B qp \ = \B q e- u -/ kT - B p e- u ^ kT \ » 0. This 
occurs in two cases: 

(1) . B q e- u «' kT » B p e- u ?/ kT , then B qp « B q e~ u ^ kT ; 

(2) . B q e- u "/ kT « B p e- u "/ kT , then B qp 
-B p e~ u ^ kT . 

Here we shall calculate the solutions only using the 
linear activation energy^ 



U q {B q ) = 17 90 (1 - Bg/B. 



q0) 



(10) 



The reason is that the calculation is very difficult using 
other activation energies. Later we shall further show 
that Eq.(10) can be used to describe the vortex motion 



in surface pinning region. 



2. B qp w -B p e~ u ' /kT 
Substituting B qp sa —B p e~ Up / kT into Eq.fjl]), we have 
^ = -e q e' u ^ kT + e qp B p e^ kT . (16) 



Substituting Eq.(10) into Eq.(16), we have 



dy 
dt 



(e q p~P q Bpe- u «l kT )y=± 



where 



(17) 



(18) 



and the initial condition B q (0) = B q o is now replaced by 
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2/(0) = 1- 

Rewriting Eq.Q as e~ u p/ kT = -(dB p /dt)/e p and sub- 
stituting it into Eq.(17l, we have 



dy 
dt 



^7dt {B *>> 



y 



l 



(19) 



general solution is 



Eq.(|19|) is a first order linear differential equation. Its 

(20) 



y(t)=P-\t) 



— [ P(t) dt + C 2 
T i J 



where 



P(t) = exp 



(21) 



and the C2 in Eq.(20l is an integrating constant that is 
determined by the initial condition y(0) = 1. The B p (t) 



in Eq.(21 1 is the internal field in the p-th pinning region. 
It has two options, i.e., Eq.(|8]) or Eq.Q. 

Finally, the internal field in the g-th pinning region can 
be obtained using Eq.(18), that is, 



S g (*)=i? 9 o{l-^m[y(i)]}. 



(22) 



In the absence of flux exchange, the parameter e qp = 
(see Eq.Q). Therefore, Eq.(22) reduces to a logarith- 
mic function B q {t) = B q0 [l - ^ x In (1 + t/r q )] . This is 
consistent with our assumption that the vortex activa- 
tion energy in the g-th pinning region obeys the linear 
law U q {B q ) = [7,0(1 - B q /B qQ ). 



V. DISCUSSION 

In this section, we shall prove that the system (Eq. ^ , 
Eq.Q) is equivalent to Maxwell's equations. We shall 
also discuss the possible flux exchange between the bulk 
pinning region and surface pinning region of a type-II 
superconductor. 



A. Equivalent with Maxwell's equation 

Because the vortices in a type-II superconductor arc 
the quanta of magnetic field, their motions must obey 
the Maxwell's equation^!, or continuity equatioiP. This 
equation gives 



d t B = ~d x (Bv) . 



(23) 



where B is average internal field (vortex density) and v 
is the vortex velocity (see Eq.Q). Let us now prove that 
Eq.@, Eq.Q are equivalent to Eq.([23|. 



As mentioned before, the vortex hopping motions in 
the smaller regions are described by different equations. 
This indicates that the right side of Eq.(23) cannot be 



directly used to describe the flux exchange between two 
contacting pinning region. To solve the problem, we may 
refer to the definition of one-sided derivatives in math- 
ematics. Let us now consider the left derivative of Bv, 
that is, 



d- (Bv) = lim 

/i->0- 



B(x + h)v(x + h)- B(x)v(x) 



(24) 



Let x + h be a point in the p-th pinning region and x 
be a point in the g-th pinning region. Therefore, Eq.(23) 
can be rewritten as 



dB p 
dt 



exchange 



-C p (B p v p - B q v q ) . (25) 



where C p is an unknown constant. 

On the other hand, the dissipation of the vortex motion 
(flux relaxation) should also cause a reduction to B p , the 
average internal field in the p-th pinning region. This 
process is usually described by an approximatiorP2 of 
Eq.|23l, that is, 



dB, 
dt 



-e„e 



-U p /kT 



relaxation 



Combing Eq.((25|) and Eq.|26|), we have 

C p SVQB pq . 



dB v 



dt 



-e„e 



-U p /kT 



(26) 



(27) 



One can see that Eq.Q and Eq.(27) are the same 
equation (Eq.([l| is used), where C p = ^^1 iL pq j A p . This 
shows that Eq.(2|) is equivalent to the Maxwell's equa- 
tion Eq.|23]). Similarly, we can also prove that Eq.Q is 
equivalent to Eq.(23). 



B. 



Flux exchange between bulk pinning region and 
surface pinning region 



In the superconductor's bu lk, the vortices are sub- 
jected to weak random pinning^* 10 * 35 ! and have a higher 
energy^. In a thin layer close to the superconductor's sur- 
face, h owever, the vortices are subjected to strong surface 
pinning^EHMl anc [ have a lower energyP Therefore, the 
bulk pinning region and surface pinning region should 
exchange magnetic flux (see Fig. 3). 

Let us now consider a cylindrical superconductor with 
a radius R. An external magnetic field is applied par- 
allel to the axis of the superconductor. Assign the bulk 
pinning region as the p-th pinning region and surface 
pinning region as the g-th pinning region. Let 6 be 
the thickness of the surface pinning region (5 is of the 
order of the superconductor's penetration depth, and 
S << R). Therefore, the area of the bulk pinning region 
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FIG. 3: (Color online) Schematic diagram of magnetic flux ex- 
change between the bulk pinning region and surface pinning 
region of a type-II superconductor. The black disks represent 
vortices. The arrows represent the vortices' hopping direc- 
tions. A p is the cross area perpendicular to the magnetic 
field in the bulk pinning region. A q is the cross area perpen- 
dicular to the magnetic field in the surface pinning region. 5 
is the thickness of the surface pinning region. 



other. Moreover, the surface pinning force is stronger 
than the interacting force between vortices because the 
surface pinning is a strong pinning. Thus, we can safely 
ignore the higher order terms in U q (B q ) and assume that 
it obeys the linear law Eq.(lO). 

3. Total internal field. — In an experiment, a measured 
internal field is usually the total average internal field 
over the entire superconductor. This field can be defined 
as 



B t = 



*P + *q 

A 



= (1 - a)B p 



aB, 



(28) 



where A = A p + A q = irR 2 is the total area of the surface 
perpendicular to the magnetic field. B p = <& p /A p and 
B q = & q /A q are the average internal fields in the bulk 
pinning region and surface pinning region, respectively. 
The geometrical factor (area ratio) 



An ( S 

q = 1 - 1 - 



.4 



R 



is A p = tt(R— 8) 2 , the area of the surface pinning region 
is A q = ttR 2 — ir(R — S) 2 , and the length of the line on 
the interface perpendicular to field is now the perimeter 
of the bulk pinning region, i.e., L pq = 2tt(R — S). 

Let R — > co, we have L pq /A p — 2/(R — 5) — > 0, and 
consequently, e pq 0. This shows that the first approx- 
imation e pq B pq w is satisfied. Thus, Eq.([7| can de- 
scribe the internal filed in the bulk pinning regio n. On 
the other hand, the bulk pinning is a we ak p inning 9 10 35 
and surface pinning is a strong pinning^IlM^ we have 
U p « U q , B p e~ u p /kT » B q e~ u "/ kT , and consequently, 
\B qp \ « \B p e- u "/ kT \ » 0. This shows that the sec- 
ond approximation |e gp i? g)9 | >> is also satisfied. Thus, 



is the weight of the contribution from the magnetic flux in 
the surface pinning region. Let R — > co, we have a —¥ 0. 



Eq.(16l can describe the internal field in the surface pin- 



ning region B q (t). 

2. Activation energy. — In the bulk pinning region, 
there is a number of options for the activation energy 
U p (B p ). However, the exact solutions can be found only 
with the linear activation energy and logarithmic activa- 
tion energypS 

In the surface pinning region, let us now prove that 
the linear activation energy (see Eq.(lOl) is more accu- 
rate than others. Recall that the vortex deformation and 
the nonlinear interaction between vortices result in the 
higher order (nonlinear) terms in the activation energy!^ 
However, the vortex deformation in the surface pinning 
region is small because the vortices are parallel to each 



Thus, Eq.(28) reduces to B t = B p . This means that, in a 



large superconductor, the contribution to the total aver- 
age internal field B t from the magnetic flux in the surface 
pinning region aB q is ignorablc. This property is com- 
pletely determined by geometry and is nothing related to 
physics. 



VI. CONCLUSION 

In a type-II superconductor, the different pinning 
regions exchange magnetic flux (vortices). This phe- 
nomenon causes the vortices in a weak pinning region 
to migrate to a strong pinning region. Thus, it influences 
the flux relaxation behavior in each pinning region. The 
flux exchange phenomenon also occurs between the bulk 
pinning region and surface pinning region of a type-II 
superconductor. But the vortices hop from the surface 
pinning region to the bulk pinning region has little in- 
fluence on the average internal field in the bulk pinning 
region. In the calculation of total average internal field, 
the contribution from the surface pinning region is sig- 
nificant in a small superconductor, but is ignorable in a 
large superconductor. 
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